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Abstract 

We review applicability of QCD factorization theorem to multiple scattering 
in deeply inelastic lepton-nucleus scattering. We show why ^4 1//3 -type nuclear 
enhancement can be calculated consistently in perturbative QCD. We derive 
the transverse momentum broadening of the leading pions in deeply inelastic 
lepton-nucleus collisions. We argue that the measurement of such transverse 
momentum broadening can provide direct information on the multiple-parton 
correlation functions inside a nucleus. We also estimate the numerical values 
of the broadening at different values of xb and Q 2 . 
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I. INTRODUCTION 



Multiparton correlation functions inside a large nucleus are extremely important and 
useful for understanding nuclear dependence in relativistic heavy ion collisions. When a 
quark or a gluon passes through the nuclear matter, it loses energy via radiation, and it 
picks up extra transverse momentum because of multiple scattering The multiple 

scattering is directly associated with multiparton correlation functions ||. Inside a large 
nucleus, multiple scattering can take place within one nucleon or between different nucleons. 
At high energy, a single hard scattering is always localized within one nucleon, and it cannot 
generate any large dependence on the nuclear size. Similarly, multiple scattering within one 
nucleon cannot provide much dependence on the nuclear size either. Therefore, the large 
dependence on the nuclear size is an unique signal of multiple scattering between nucleons. 
Measurements of such anomalous dependence on the nuclear size for any physical observable 
will provide direct information on multiparton correlation functions in a nucleus. 

However, there are many multiparton correlation functions in QCD ||. It is therefore 
important to identify physical observables which depend only on a small number of mul- 
tiparton correlation functions. Otherwise, data will not be able to separate contributions 
from different correlation functions. It was pointed out recently in Ref. || that the A 1 ^ 3 - 
type enhancement of the transverse momentum broadening of Drell-Yan pairs and the jet 
broadening in lepton-nucleus deeply inelastic scattering (DIS) are good observables, because 
at the leading order of a s , these two observables depend only on one type of multiparton 
correlation functions. It represents the correlations between hard quarks and soft-gluons, 
and has the following operator definition |)]||, 

T^^ 2 ) = I ^e***y-y&*0te -y-)0{y- 2 ) 

x (Pa\M0) £ F\{y-) F +a (y~) ^(y") \P A ) , (1) 

with A the atomic weight, and q is the quark flavor and F is the gluon field strength. In 
Eq. flU), x is the effective momentum fraction and /i represents the scale-dependence of the 
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correlation functions. 

According to QCD factorization theorem all information on the identity of the 

target is factorized into the matrix elements. The A 1 / 3 -type enhancement must be a property 
of the relevant nuclear matrix elements, like the one in Eq. (|T]). Experimentally, the A 1 / 3 - 
type nuclear enhancement have been observed in the transverse momentum broadening 
of Drell-Yan pairs in hadron- nucleus collisions |§. Theoretically, Eq. ([!]) shows how such 
enhancement can occur, through integrals over the relative positions of the fields, i/i, that 
appear in the expectation values. For example, enhancement can occur when the two quark 
fields are close together and the two gluon fields are close together, and they pair off into 
color singlets that are separated by no more than a nucleon diameter in the i/i integrals; 
but the two pairs are separated by a distance that varies up to the nuclear size. In this 
manner, aside from an overall factor of A which reflects growth with the nuclear volume, we 
anticipate an A 1 / 3 nuclear enhancement fL0| . 

Direct measurement of the multiparton correlation functions is of great importance for 
testing the generalized QCD factorization theorem ||, which allows us to apply QCD per- 
turbation theory for studying the collisions involving nuclei. It was observed in Ref. |5| 
that measuring the jet transverse momentum broadening in DIS, A(p^), provides a direct 
measurement of the quark-gluon correlation functions, TA(xb, Q 2 ), inside a nucleus. At the 
lowest order, the jet transverse momentum broadening in DIS can be expressed as || 

^ 2 a s (Q 2 ) E q e 2 q T q A F (x B ,Q 2 ) 

3 Z q e 2 q A (x B ,Qi) ' 1 } 

where Pt is the transverse component of the jet momentum in the photon- nucleus center 
of mass frame in DIS. In Eq. (|2|), x B is the Bjorken variable and Q 2 = —q 2 with q 2 is the 
invariant mass of the virtual photon in DIS. 

However, other than a future HERA with a nuclear beam, existing fixed target facilities 
can not provide good measurements on jets in lepton-nucleus DIS. It is the purpose of this 
paper to show that by measuring the leading pions in DIS and their transverse momentum 
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broadening, fixed target facilities, such as CEBAF [[11]], can provide good information on 
the quark-gluon correlation functions defined in Eq. (JT]). 

This paper is organized as follows. In the next section, we derive the transverse momen- 
tum broadening of the leading pions in lepton-nucleus DIS. We demonstrate analytically 
the direct correspondence between the transverse momentum broadening of the leading pi- 
ons and multiparton correlation functions. In Sec. Ill, we estimate the numerical values of 
such broadening at different values of xb and Q 2 , by using the model of the quark-gluon 



correlation functions introduced in Ref. |10| . We also derive explicit relations of the trans- 
verse momentum broadening for 7T ± and 7r°, and discuss the xb and Q 2 dependence of the 
correlation functions. 

II. TRANSVERSE MOMENTUM BROADENING OF PIONS IN DIS 

Consider leading pion production in the deeply inelastic lepton-nucleus scattering, 

e(fci) + A(p) ^ e(k 2 ) + tt(£) + X , (3) 

where ki and k 2 are the four momenta of the incoming and the outgoing leptons respectively, 
p is the momentum per nucleon for the nucleus with the atomic number A, and i is the 
observed pion momentum. In order to extract useful information on multiple scattering from 
the pion production, it is natural to think that the A-dependence of the differential cross 
section da e A^ e TrX /dxsdQ 2 d£ 2 n is the physical observable to study, because the transverse 
momentum spectrum is more sensitive to multiple scattering ||. But, other than a future 
HERA with a nuclear beam, existing fixed target facilities can not produce good data on the 
production of leading pions at large transverse momenta in DIS. On the other hand, at the 
small transverse momentum l\, the differential cross section da j dxBdQ 2 dl\ is proportional 
to 1/i 2 ^, and is not perturbatively stable. A nontrivial all order resummation of the large 
logarithms log(Q 2 /£^) is necessary in order to get a reliable prediction of the transverse 
momentum spectrum at small Gp [I2|. However, the inclusive moments of the transverse 
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momentum spectrum / d£\{£\) n da e A^enX / ' dx sdQ 2 'di\ with n > 1 are perturbatively stable, 
and also enhance the information in the region of large l\ y where the multiple scattering is 
most relevant. Therefore, the A-dependence for the moments of the transverse momentum 
spectrum of the leading pions is a good observable for extracting information on multiple 
scattering. 

We define the first moment of the transverse momentum spectrum - the averaged pion 
transverse momentum squared, as 



<4) eX = <4=S/=^f; (4) 



, 2 \i . I _ / , j , 2 , 2 d<J e A^eirX J do e A^eX 

dxBdQ 2 d£^ J dxsdQ 2 

where da e A^ e x / dxsdQ 2 is the total inclusive DIS cross section normalized by the atomic 
weight A. In Eq. (f|), the Bjorken variable Xb = Q 2 /(2p • q), and q = k\ — k 2 is the four- 
momentum of the virtual photon. The transverse momentum of the leading pion depends 
on our choice of the frame. In this paper, we choose the photon-nucleus frame, and choose 
the z-axis along the momentum line of the nucleus and the virtual photon. 

To separate the multiple scattering contribution from the single scattering contribution, 
we define the nuclear broadening of the transverse momentum square as 

a(4) = (f T r A - (f T y N . (5) 

As we will demonstrate below, the nuclear broadening of the transverse momentum square 
defined in Eq. @ can be parameterized as |ll3| , |T4| 



A(4) = a + bA 1 / 3 , (6) 

where bA 1 ^ 3 represents the contribution directly from the multiple scattering which is ex- 
plicitly proportional to the nuclear size (oc A 1 / 3 ) and term a includes all contributions from 
the localized hard scattering as well as those suppressed by high power of 1/Q 2 - In principle, 
the parameter a in Eq. (0) can also depend on the atomic weight A. But, as we will explain 
below, its dependence on A should be very weak (e.g., A a dependence with a ~ ±0.02). In 
this paper, we are interested in the second term in Eq. @, and we show that experimental 



measurement of the constant b in Eq. (|6|) will provide direct information on the quark-gluon 
correlation functions shown in Eq. ([[]). 

To derive explicit expressions for the a and b in Eq. (|^), we expand both numerator and 
the denominator in Eq. in terms of a power series of 1/Q 2 . In addition to the terms 
of leading power, we keep only the power suppressed terms that are explicitly enhanced by 
a factor A 1 / 3 . Since the denominator, da e A^ex/ dxsdQ 2 , is the total inclusive DIS cross 
section, the operator product expansion (OPE) allows us to expand it in terms of the power 
series of 1/Q 2 [||,@, 



dxsdQ 2 



j dx a /A(x,/i 2 



da. 



(0) 

ea— *eX 

dxsdQ 2 



(x B /x,jj 2 /Q 2 ,a s (iJ 2 )^ 



= D 



(0) 



1 + 



,Q 2 



(?) 



where /i is the factorization and/or renormalization scale, 4> a /A{%, Z^ 2 ) is the leading twist 
parton distribution of flavor a in a nucleus normalized by the atomic weight A, and 
dcr^^ex/ 'dxsdQ 2 is a perturbatively calculable short-distance hard part, which is indepen- 
dent of the nuclear medium. The total DIS cross section normalized by the atomic weight 
A is an inclusive quantity and depends only on one hard scale Q 2 . The QCD factorization 
theorem allows us to factorize each term in Eq. (|7|) into a convolution of a short- 

distance coefficient function and a corresponding matrix element. The 1/Q 2 term in Eq. (^) 



can be expressed in terms of twist-4 matrix elements fL6 |. As demonstrated in Refs. |p|,[17|, 
integration of all position variables (yt) of the field operators, which define these twist-4 
matrix elements, are bounded by 1/xbP due to the oscillating exponential. Let m and r be 
nucleon mass and radius, respectively. If 1/xbP < 2r(m/p) (i.e. xb > l/2mr ~ 0.1), all the 
fields in these twist-4 matrix elements are bounded within the size of individual nucleon. 
Therefore, the power suppressed terms in Eq. ([TJ) is of the order of 0(1/Q 2 ) if xb > 0.1, not 
OiA^/Q 2 ). 

On the other hand, the numerator in Eq. ([|) depends on two physically measured 
hadrons: the pion and the nucleus, and therefore, the OPE alone cannot be used to ex- 
pand the numerator. However, the generalized factorization theorem introduced in Ref. f| 
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can be used to factorize the numerator up to the 1/Q 2 power corrections. Similarly to the 
Eq. (1) of Ref. flQ|l , we expand the numerator in Eq. (|J) as 

/ dfr 4 dx B dQHl 2 T = ? ® C ^ ecX ( Xb/x > Z = i/p " Q2 ) ® Dc ^ (0) 

a,c 

+<Pa/A{x) ® <?£L eC X (W^, Z = Q 2 ) ® 4-tt(*) 

+ ... 

ssltf>+ltf > + #£> + ..., (8) 

where ... represents terms further suppressed in 1/Q 2 , ® represents the convolution over 
partons' momentum fractions, and explicit dependence on the factorization and/or renor- 
malization scale is suppressed for simplicity. In Eq. (|8|), C^°\ C^ 2 \ and are perturba- 
tively calculable short-distance hard parts ||, and D c ^ and d c ^ are twist- 2 and twist-4 
parton-to-pion fragmentation functions, respectively. <p a /A are nuclear parton distributions, 
and the T a /A in Eq. (||) represents the twist-4 quark-gluon correlation functions (e.g., the 
one defined in Eq. ([I])). Both 4>a/A an d ^ are normalized by the atomic weight A. 

Effective nuclear parton distributions <f> a /^(x) have been measured [IB-CT] and are known 



to have the nuclear shadowing for small x (x ~ 0.1), the EMC effect for intermediate x 
values (0.3 ~ x ~ 0.7), and Fermi motion effect for the large x region. By fitting recent high 



precision DIS and Drell-Yan data on various nuclear targets, Eskola et al. |22| extracted 
a set of effective scale dependent nuclear parton distributions. In order to estimate the 
A-dependence of the nuclear parton distributions, we introduce a parameter a(A, i B ) as 



In 

F 2 A (x B , Q 2 ) = A a ^ F 2 N (x B , Q 2 ) , or a(A, x B ) = - 



F A (x B ,Q 2 ) F 2 N (x B ,Q 



In [A] 

(9) 

where F A (x B ,Q 2 ) is a nuclear structure function normalized by A and F 2 N (x B ,Q 2 ) is an 
isoscalar nucleon structure function. Using the lowest order expression for the structure 
function: F 2 (x B ,Q 2 ) = J2 q e l x b (pg{x B , Q 2 ) with e q being a quark's fractional charge and (f) q 



the effective nuclear parton distributions of Ref. [22]], we plot the parameter a(A, xb) as a 



function of xb for A=12, 64, and 207 at Q 2 = 4 and 9 GeV 2 in Fig. |T]a and |I]b, respectively. 
Fig. ^ shows that although the exact xe-dependence and the A-dependence of the structure 
functions are non-trivial, the overall A-dependence of the structure functions is limited to 
A a with a value of a ~ ±0.02 for xb values between two vertical lines in Fig. [IJ, which 
are relevant to pion production in this paper. Actually, as shown in Fig. [j], the absolute 
value of a is limited to 0.05 for any practical value of xb accessible at fixed target energies. 
Such A-dependence in the structure function as well as in the 4> a /A{x) is much weaker than 
the v4 1//3 -dependence caused by multiple scattering. Since the parton-to-pion fragmentation 
functions, D c ^ n and <i c ^7r, should not have explicit A 1 / 3 -dependence, the term is of the 
order of 0(1/Q 2 ), not 0(A 1 ^ 3 /Q 2 ). On the other hand, the term depends on the multi- 
parton correlation function in a nucleus T a /A(x), and therefore, it will have an A 1 / 3 -type 
enhancement |10| . 



Substituting Eqs. ([j]) and ([8]) into Eq. (|), and keeping terms up to 0{A 1 ^ j 'Q 2 ), we 
obtain 



M 0) W 2 , 



A +0[- 2 ) • (10) 



Similarly, for a nucleon target, we have 

<4r*§ + o(|). (id 

Substituting above Eqs. ( |To|) and (|ll|) into our definition of the nuclear broadening of the 
transverse momentum square in Eq. (|5|), we derive 

Hf „ I A" 

,<»> 

where the power suppressed terms 0(A°/Q 2 ) should be very small due to the cancellation 
between {&p) eA and {£ 2 ^) eN and the fact that they are not enhanced by the A 1 / 3 . The first 
term in Eq. flT2"D would be exactly equal to zero if the normalized effective nuclear parton 
distributions 4> a /A( x ) = 0a/Ar(^)- However, due to the well-known nuclear effects in the 
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1 



tr(0) rr(0) 
A __ H N 

D (o) D (o) 



+ W + {^' (12) 



parton distributions, the first term in Eq. (|12D does not have to vanish. But, because of very 
weak A-dependence of nuclear parton distributions shown in Fig. [I], we identify this term 
with the a in Eq. (IS), 



rj-(O) rr(0) 

n A H N 



D (o) D (o) • 

U A N 



(13) 



It is the second term in Eq. (|12D that is responsible for the v4 1//3 -type dependence of the 
nuclear broadening. 

We introduce A(£j)i/3 as 

A(4)i/ 3 = 6^ « ^ 



t / , (0) 

«2 /)2 aCr eA^e7rX / aCr eA^eX 



(14) 



dx B dQ 2 dl% I dx B dQ 2 

where superscript "(D)" represents the double scattering contribution, and "(0)" stands for 
the leading power contribution. Our formalism in Eqs. (|13|) and (T3[) can be verified from 



the A-dependence of the measured transverse momentum broadening. 

From the factorized formula in Eq. (H), the double scattering contribution depends 
only on the twist-2 parton-to-pion fragmentation function D(z). Therefore, the differential 



cross section in Eq. ([14 ) for the pion production in DIS can be expressed as 



da {D) r da {D) dr 



dx B dQ 2 di\ ? J dx B dQ 2 dp 2 CT ' Dc ^ z) ' z i ' (15) 
where J2 C runs over all parton flavors c, z = £/p c is the momentum fraction carried by the 
produced pion, and D c ^(z) is the fragmentation function for the parton of momentum p c to 
become a pion of momentum i. In Eq. (|I5D, da^_^ ecX / dx B dQ 2 dp 2 T is the double scattering 
contributions to the differential cross section to produce a parton of momentum p c in DIS, 
and it can be written as 



a<J eA^ecX 



d Pl = ^^^^{k x M)W^\x B M\Vc) , (16) 
1 Sir x B s z Q z p 



dx B dQ 2 dp 2 CT 

where s = (p + ki) 2 is the total invariant mass of the lepton-nucleon system. In Eq. ( |T5D 
the leptonic tensor L^ v is given by the diagram in Fig. [|a, 
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L^(k 1 ,k 2 ) = -Tr( 1 -k 1 Yl-hY) 



(17) 



and WS™ is the hadronic tensor due to double scattering, which is given by the diagram 
shown in Fig. |2|b. For comparison, the lowest order differential cross section for producing 
a parton of momentum p c due to single scattering is given by 



da. 



(0) 

eA^ecX 



dxsdQ 2 dp 2 



dp. 



CT 



da 



(0) 



dp, 



(18) 



dxBdQ 2 

where da 1 ^ j dxsdQ 2 is the leading order inclusive DIS cross section, which appears in the 
definition of A(4) 1/3 in Eq. (0). 

The complete double scattering contributions to the production of a leading quark of 
momentum p c in DIS at the leading order of a s come from the diagrams shown in Fig. |||, 
which represent the final-state interactions, plus the same order diagrams involving initial- 
state interactions shown in Fig. |j. In the following, we first discuss the role of initial-state 
interactions, and argue that although required by gauge invariance, they contribute to the 
a term of Eq. (|6]) only. 

After the collinear expansion of the parton momenta, the gluon interactions on the initial 
quark lines, as shown in Fig. [|, can be reduced into two categories: the long-distance and 
the short-distance contributions due to the fact that every "on-shell" propagator can have 
a pole contribution as well as a contact contribution |T7|]. For example, a quark propagator 
of momentum k can always be written as 

ij ■ k ry • k 



Z7 • n 

+ 777 , 



(19) 



k 2 k 2 2k ■ n 

where k 2 = and n M is any auxiliary vector with k - n ^ 0. The first term in the right-hand- 
side of Eq. ([I9| ) corresponds the pole contribution, while the second term is the contact 
contribution |T/]. Attaching one gluon to the initial quark line introduces a quark prop- 
agator, and this propagator will have both the pole and contact contributions. The pole 
contribution is long-distance in nature, representing the interactions between the quark and 
the gluon long before the hard collision between the quark and the virtual photon. The 
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pole contribution is partially responsible for the A-dependence of the leading-twist parton 
distributions in a nucleus and therefore, is one of the sources for the weak A-dependence 
appeared in the a term in Eq. (^|). 

Because of the nature of the contact term, its contribution is localized in space, and 
therefore, it does not contribute to the A 1 / 3 -type of nuclear enhancement ||. Since the short 
distance contributions of Feynman diagrams in Fig. § have at least one propagator given by 
the contact term, these diagrams do not contribute to the A^-type of nuclear enhancement. 
At the same time, the contact term of the initial-state interaction is important for the gauge 
invariance of the complete double scattering (twist-4) process |[L7|| , but only at order A . 

Although the Feynman diagrams with the final-state interactions in Fig. |3] are not gauge 
invariant in general, their contributions to the A 1 / 3 -behavior of nuclear enhancement are 
observable and hence gauge invariant ||10|| . The three lowest order diagrams in Fig. |3] are all 



convoluted with the same two-quark-two-gluon matrix element through three independent 
parton momenta. The leading power contributions of these diagrams come from the region 
phase space when all partons are collinear to the direction of the target, as shown in Fig. |3|. 
In order to convert the gluon field A + to the corresponding field strength F +± in covariant 
gauge, we expand the gluon momenta in an extra transverse component kx in Fig. [3] [IU|. All 



three Feynman diagrams in Fig. ^| have two potential poles due to the interactions between 
the final-state quark and gluons. These two poles (which are not pinched) can be used to 
carry out the momentum fraction dx\ and dx2 integrals P,|10||. After taking the poles, the 
hard-scattering factor (the interaction between the virtual photon and the initial quark) 
and the final-state interaction between the "on-shell" outgoing quark and the gluon are 
separately gauge invariant. 

To derive the leading contribution from the double scattering, we follow the derivation 
in Ref. ||. By taking the poles to fix the integrations of momentum fractions X\ and X2, we 
derive contributions of all three diagrams in Fig. to the hadronic tensor Wffl in Eq. ([16]) 
as 
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Pc 



d 2 k elkT -y 1T J p+( Vl -y 2 ) 



J_ r dy dy l dy 2 d 2 y T 
4vr J 4tt (2tt) 2 (2tt) 2 

x (p A \$Ml + A + (y^,y 1T ) A + (y 2 )^ q (y-) \p A ) 



x(2n)6(y^-y-)(2n)e(y^e i ^ + y- 
2tt 

W~q) 3 """ 



x 



i 



7T 



dp; 



2 



(20) 
(21) 



_(2p-g)3 

In Eq. (pi|), T^(:rB, /i 2 ) is denned by Eq. (p with the factorization scale /x 2 . Hff'(p c ) 
represents the spinor trace of the partonic part of the double scattering and is given by 

H\H\pc) = -7r 2 a s a em e 2 q Tr [7 • p7 M 7 • (xp + q)Yl • Pd P l • (xp + q)^f v \ p P p a , (22) 

where a color factor 1/2N C with N c = 3 was included. Using the definition z = £/p c , see 
Eq. fl25|) below, we can reexpress the momentum p c in S'{p 2 ) in terms of the momentum i 
of the observed pion, 

d 



C>2 



d{% K" z2 ' 



z 4 8\e T ) . 



(23) 



Combining Eqs. (|15"D , (|16D, (pT|) , and (p3j), we obtain 



<i<r. 



7T 



dx B dQ 2 di\ 8nx 2 B s 2 Q 2 (2p ■ q) 3 



E / ^ Dq^^z) T^ F (x B , Q 2 ) L^ v (£ / z) [-£'(4) 



(24) 



where the runs over all quark and antiquark flavors, and z = £/p q is the momentum 
fraction of the quark flavor q carried by the observed pion. At the lowest order in a s , hard 
gluon initiated double scattering, which is proportional to the four-gluon correlation function 
T gF , does not contribute HJlO[ . 

In the photon-nucleus frame, we choose the target momentum p along the — iT-axis, such 
that p^ = (p° , p x , p y , p z ) = (P, 0,0, — P), and only p~ = (p° — p z )/\/2 is nonvanishing. In 
this frame, we have the following expression for the momentum fraction z = £/p c , 

t _ e + p-t 

Pc pi 
p-t 



z = 



P'Pc 



p-q 



(25) 
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In deriving the last equation, we used p c = xbP + q and p 2 ~ 0. 

After working out the algebra for L^Hffi, and substituting Eq. ( p^ ) into Eq. (fH|), we 



obtain 



A(4) 1/3 (4nn) 



3 E q e 2 q q A (x Bj n 2 ) 



(26) 



where D g 



-7T 



(z) are the quark-to-pion fragmentation functions, and q A 



(X B)/ U 2 ) = </> q / A (x B , P?) 



are the effective quark distributions in the nucleus normalized by the atomic weight A. In 
Eq. fl26D , z m i n = p ■ £ m in/(p ■ q) is defined in terms of £ m i n , which is the momentum cut for 
the pions measured in the experiment. One can choose £ m j n to be large enough to ensure 
that the observed pions are from the fragmentation of energetic quarks. 

In Eq. (p6|), the jj 2 in a s {ji 2 ) is the renormalization scale, and the fi 2 in T a f {xb, /i 2 ) and 
q A (xB, A* 2 ) are the factorization scale. Since the short- distance part of the transverse momen- 
tum broadening defined in this paper is an inclusive and infrared safe calculable quantity, 
the scale /i 2 should be chosen to be the order of the physically measured momentum scales. 
As shown in Eq. (pE|), the only physically measured momentum scales for the transverse 
momentum broadening are Q and £ m i n . Since \£ m i n \ are of the same order as y/Q 2 ', we choose 
/i 2 = Q 2 for the numerical calculations below, and we believe that such a choice will not 
result into the large logarithmic high order corrections. 

Since the quark-to-pion fragmentation functions have been measured ||24|| , Eq. fl26|) shows 
that the transverse momentum broadening of pions provides direct information on the quark- 
gluon correlation functions inside a nucleus. The measured size of the transverse momentum 
broadening, A(£ 2 -) 1 / 3 , depends on the choice of the momentum cut £ min (or equivalently z m i n ). 
Our predictions should be more reliable for the leading pions, or pions with relatively large 
momenta. By measuring the transverse momentum broadening for tt^ 1 and tt°, and keeping a 
reasonable large value of z m { n , we can extract the quark flavor dependence of the correlation 
functions. 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



The numerical values of the pion transverse momentum broadening, A(£y) 1 / 3 (£ min ), de- 
pend on the explicit functional form of the quark-to-pion fragmentation functions and the 
quark-gluon correlation functions. In this section, without assuming any specific form for the 
quark-gluon correlation functions, we derive relations for the transverse momentum broad- 
ening between 7r =1= and ir° based on isospin symmetry and the charge conjugation invariance. 
Furthermore, by using the simple model proposed in Ref. Q for the quark-gluon correlation 
functions, we explore both the normalization and functional dependence of the momentum 
broadening. 

Although new parameterizations of the quark-to-pion fragmentation functions were ob- 



tained recently f24 |, we will use the simple parameterizations of Ref. |£5] for the following 
analytical derivations and discussions on the general features of the momentum broadening. 
Later, when we present our figures for the numerical values of the momentum broadening, we 



will use the parameterizations from both Refs. |25] and p6fl , and demonstrate the similarities 
and differences. 

Like the parton distributions, the quark-to-pion fragmentation functions have scaling 
violation (or Q 2 -dependence). To simplify our discussion on the general features of transverse 
momentum broadening, we ignore the scaling violation of the fragmentation functions, and 



take for the fragmentation functions the input distributions of Ref. [251, which are given as: 



D+{z) = 




U d — 


U 


= Df = 


D-(z) = 


d: = 


u d — 


D- + 

u 


= Df = 


D\z) 


= Df - 


- Df = 


'Df 


= Df = 



1-z 2 



:i 



z 



Az 
1 - z 



Az 



(27a) 
(27b) 
(27c) 



[D + (z) + D (z)}/2. For the strange quark fragmentation functions, we use 
DZ^ = = = D~ . Notice that the fragmentation functions given in Eq. (P?D, 



with D"(z 

T + 



as well as those given in Ref. p4| , p6| , violate Gribov-Lipatov reciprocity, which requires 
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that the power of (1 — z) must be even. However, since it is not our purpose to invent 
better fragmentation functions in this paper, we will first use the fragmentation functions 
in Eq. ( p7|) in our numerical calculation to illustrate the size and the general features of the 
transverse momentum broadening. Our predictions can be easily adjusted for other sets of 
fragmentation functions. 

In order to evaluate the transverse momentum broadening, Eq. fl26|) requires the moments 
of the fragmentation functions. We therefore introduce 

D^\n,z min )= f dzz n D^(z) , (28) 

^-min 

where i = +,— and for 7r + , 7r~ and 7r°, respectively. Using the fragmentation functions 
defined in Eq. fl27|), we obtain the following identity, 

D + (n, z min ) - D°(n, z min ) = D°(n, z min ) - D~(n, z min ) , (29) 

for all n. This identity is useful for the following discussions on the general features of the 
transverse momentum broadening. From Eq. fl2"6|), the transverse momentum broadening, 
A(^T)i/3(4mn), depends on the second moments D^(2, z m - m ), which are plotted in Fig. |5|. 

Given the quark-to-pion fragmentation functions, we derive the following general rela- 
tions between the transverse momentum broadening of 7r + , 7r~, and tt° particles: 

A(4)vs - A<4>& = A(4)i/° 3 - A<4)fc , (30a) 

and 

a(4)-; 3 > A(f T )f /3 > a(4k; 3 , (sob) 

independent of the details of the correlation functions T^ F . Eq. ( |30a| ) is a result of the 
isospin symmetry of the fragmentation functions given in Eq. ( |27D and the identity given 
in Eq. (p9|). The inequality given in Eq. ( |30b| ) results from the fact that the contribution 
of each quark flavor is weighted by the square of the quark's fractional charge, e^, and the 
relation D+(z) > D°(z) > D~(z). 
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Since the correlation functions T A F are not known, in order to obtain numerical estimates 
of the transverse momentum broadening of pions, we adopt the following model for quark- 
gluon correlation functions flJTOfl: 

T A (x,Q 2 ) = X 2 A^q A (x,Q 2 ) , (31) 

where q A (x) is the effective twist-2 quark distribution of a nucleus normalized by the atomic 
weight A, and A is a free parameter to be fixed by experimental data. This model was 
proposed after comparing the operator definitions of the four-parton correlation functions 
and the definitions of the normal twist-2 parton distributions. Substituting Eq. ([31]) into 
Eq. (p6|) , we have 

2 47r 2 a s {Q 2 ) 1/3 ^ J2 q e 2 q g A (x B ,Q 2 )D q ^(2,z min ) 

A(Mi/3= ^^ A A — T^H^m — ■ { ) 

From Eq. ([32]), we can deduce a very simple formula for the transverse momentum 
broadening of 7T° by using the fragmentation functions given in Eq. ( p7cj ): 

A(4)v3 « 4?r ^3 (g2) ^ 1/3 A 2 D°(2, z min ) . (33) 

In deriving Eq. (|33"D, we neglected the strange quark contribution which is much smaller. 
This simple expression shows that A(£y)w 3 has a scaling behavior as xb varies. The ap- 
proximate xg-scaling behavior of A(£j)y 3 is a direct consequence of the model for T A F {x), 
given in Eq. ([H]). If the future experimental data on the transverse momentum broadening 
of 7T° shows a strong violation of the x^-scaling, the model for the quark-gluon correlation 
function TA,(xb), given in Eq. fl3"T]), will have to be modified. In addition, Eq. ([53]) also 
shows that Q 2 dependence of A(£^)^ 3 is mainly from the Q 2 dependence of a s (Q 2 ). Since 
the Q 2 dependence of a s (Q 2 ) is already known, we can learn the Q 2 dependence of the cor- 
relation function T a f (xb, Q 2 ) by measuring the Q 2 dependence of A(£y)^ 3 . If the measured 
Q 2 dependence of A(£^)^J 3 is consistent with the Q 2 dependence of a s (Q 2 ), it means the 
correlation function T a f (xb,Q 2 ) has a similar scale-dependence as a normal parton distri- 
bution. By examining the x^-scaling property and Q 2 dependence of A(£^)^° 3 , we can test 
whether the model for T a f {xb,Q 2 ), given in Eq. (pl|), is reasonable. 
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It should be emphasized that our general conclusions given in last paragraph, as well as 
numerical estimates given below, are the consequences of our lowest order calculations. High 
order corrections in a s and/or in inverse powers of Q 2 can certainly change the dependence 
on and Q 2 , as well as the absolute values of our numerical estimates. But, because of 
the inclusivity and the infrared safety of the transverse momentum broadening, and the lack 
of the large logarithms of the ratios of the physically measured scales, we believe that the 
change should not be very dramatic, and that our predictions should not be off by orders of 
magnitudes. 

In the following, we obtain our numerical estimates of transverse momentum broadening 
for 7r + , 7r~, and 7r° particles by evaluating Eq. (|32|). We normalize A(£|,) 1 / 3 by A 1//3 , and 
plot our results in Fig. ^ to Fig. [7]. In view of the limited range of Q 2 available at CEBAF, 
we neglect the scaling violation of the fragmentation functions for our numerical estimates. 
We used the simple fragmentation functions given in Eq. (|27|) to produce Fig. || to Fig. [7|. 
More realistic parameterizations for quark-to-pion fragmentation functions were obtained in 
Refs. [21 1 and [^(J. In particular, those in Ref. are extracted from DIS data. We present 



both results for A(4)^ /3 /A 1 / 3 by using two different sets of the fragmentation functions 
at Q 2 = 4 GeV 2 in Fig. [§] to illustrate the uncertainties in our numerical estimates due to 
different choices of the fragmentation functions. Although numerical details of the A 1 ^ 3 -type 
transverse momentum broadening depend on the fragmentation functions used, the overall 
shape and magnitude of the broadening are consistent. The most important feature is that 
the size of the A 1//3 -type broadening is large enough to be measured experimentally fL3,14 



In obtaining our numerical results, we used CTEQ3L parton distributions of Ref. [27 
for the quark distributions in the nucleons. We took an average of the quark distributions 
of the proton and the neutron for the effective quark distribution in the nucleus, q A {xB)- 
If the A-dependence of the nuclear quark distributions is independent of the flavor, from 
Eq. (0), such nuclear dependence can be factorized and canceled between the numerator 
and denominator. However, there is no obvious reason why the nuclear dependence of quark 
distributions is flavor independent. But, as we discussed above, in the range of x probed 
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in these experiments the effective nuclear parton distribution should have extremely weak 
nuclear dependence in comparison with the v4 1//3 -type enhancement that we discussed in 
this paper. In order to verify this feature, we used both q A (xs) = q N {xs) and the realis- 



tic parameterizations for q (xb) given in Ref. |22[ to calculate the transverse momentum 



broadening A(£^)^ 3 /A 1 / 3 . We found no noticeable difference for the transverse momentum 
broadening presented in Figs. H, |7|, and |], which were obtained with q A (xs) = q N (xb)- 

Figs. ||a and ||b shows A(£| I )i/ 3 /y4 1 / 3 as functions of z m - m at Q 2 = 4 GeV 2 , and 
xb = 0.2 and 0.3, respectively. Figs. [7|a and 0b show A(£^)i/ 3 / A 1 ^ 3 as functions of z min 
at Q 2 = 9 GeV 2 , and xb = 0.2 and 0.3, respectively. In these figures, the relations 
A<4)$ > A(£ 2 T )f /3 > A(4>-3 and A<4)^ 3 - A(£ 2 T )f /3 = A(£ 2 T )f /3 - A<4)^ 3 are clearly 
demonstrated. Furthermore, we have xb- scaling for A{P T )f /3 , and A(£ 2 T )f/ 3 / A(£ 2 T )lj 3 « 2 
for small z m i n . The approximate ratio, A (4) 1/3 /A (4) 1/3 ~ 2, is a result of the isospin av- 
eraged targets, and can be easily verified by using the fragmentation functions in Eq. (p7|). 
keeping only the valence quark contributions and z min <^ 1. 

Comparing Figs. [J]and[7|, it is evident that the absolute sizes of the A(4)i/3 decrease as 
Q 2 increases. This is caused by the Q 2 dependence of the a s (Q 2 ) in the overall factor of the 
transverse momentum broadening. Due to the available energy at CEBAF, we plotted the 
broadening at Q 2 = 4 GeV 2 in Fig. || Even though Q 2 = 4 GeV 2 may not be large enough 
to apply for the fragmentation analysis, our results can be tested at CEBAF with its future 
upgrade, and can be easily tested at RHIC with its future option of electron-ion collider. 

From Figs. § and 0, we noticed that the identity, A(4)^ + 3 - A(4)^° 3 = A(£%)f /3 - 
A (4) 1/3 j has some weak violation on the value of Xb, which is due to the fact that the 
number of sea quark increases at smaller xb in comparison with the valence quarks. 

The value of A 2 for the quark-gluon correlation function T^p, given in Eq. (|31~D, has 
not been well determined. It was estimated in Ref. [I0[] by using the measured nuclear 



enhancement of the momentum imbalance of two jets in photon-nucleus collisions [|T3|,p^4 
and was found to be the order of A 2 ~ 0.05 — O.lGeV 2 . However, the nuclear enhancement 
of the Drell-Yan transverse momentum [p8j^9| indicates a much smaller value of A 2 ||. 



Although the momentum imbalance of the di-jet data depends on the final-state multiple 
scattering, and the Drell-Yan data is an effect of initial-state multiple scattering, the leading 
order calculation indicates that A 2 from the two data sets to be the same ||. Therefore, 
this work on the transverse momentum broadening in pion productions should provide a 
new and independent measurement of the size of the quark-gluon correlation functions, and 
the value of the A 2 , which is extremely important for understanding the nuclear dependence 
and the event rates in the relativistic heavy ion collisions. Since both di-jet momentum 
imbalance and the transverse momentum broadening are caused by the final state multiple 
scattering, we used A 2 = 0.05 GeV 2 for Figs. ^| to [8]. Any change in A 2 will only change the 
overall normalization factor in A(£ 2 ^) 1 / 3 /A 1 / 3 , and hence does not affect the general features 
of A(4) 1/3 /A 1 /3. 

We conclude that the transverse momentum broadening of leading pions in deeply inelas- 
tic lepton-nucleus scattering is an excellent observable to probe the parton correlation func- 
tions in the nucleus. Measuring the A 1//3 -type transverse momentum broadening A(£ 2 -) 1 / 3 
in the leading pion production in DIS provides an independent test for the existing model of 
the quark-gluon correlation functions More importantly, by measuring such broadening, 
we can directly measure the x and Q 2 dependence of the correlation functions, which is very 
useful for predicting the event rates ^0] and for understanding the nuclear dependence in 
relativistic heavy ion collisions. 

ACKNOWLEDGMENT 

We thank G. Sterman for very helpful discussions. This work was supported in part by 
the U.S. Department of Energy under Contract No. DE-AC02-98CH10886, and Grant Nos. 
DE-FG02-87ER40731 and DE-FG02-96ER40989. 



19 



REFERENCES 

A.H. Mueller, in the proceeding of the Workshop on CEBAF at Higher Energies, ed. N. 
Isgur and P. Stoler (Newport News, USA, 1994). 

[2] M. Gyulassy and X.-N. Wang, Nucl. Phys. B420, 583 (1994); X.-N. Wang, M. Gyulassy 
and M. Pliimer, Phys. Rev. D51, 3436 (1995). 



[3] R. Baier, Yu.L. Dokshitzer, A.H. Mueller, and D. Schiff, |hep-ph/9907267| , and references 
therein. 

[4] G.T. Bodwin, S.J. Brodsky and G.P Lepage, Phys. Rev. D39, 3287 (1989). 

[5] M. Luo, J.-W. Qiu and G. Sterman, Phys. Lett. B279, 377 (1992); M. Luo, J.-W. Qiu 
and G. Sterman, Phys. Rev. D50, 1951 (1994). 

[6] X.-F. Guo, Phys. Rev. D 58, 114033 (1998). 

[7] J.C. Collins, D.E. Soper, and G. Sterman, in Perturbative Quantum Chromodynamics, 
edited by A.H. Mueller, 1989 (World Scientific, Singapore); and references therein. 

[8] J.-W. Qiu and G. Sterman, Nucl. Phys. B353, 105 (1991); B353, 137 (1991). 

[9] P. Bordalo et al, Phys. Lett. B 193, 373 (1987); D.M. Aide et al, Phys. Rev. Lett. 66, 
2285 (1991). 

[10] M. Luo, J.-W. Qiu and G. Sterman, Phys. Rev. D49, 4493 (1994). 

[11] S.J. Brodsky, in the proceeding of the Workshop on CEBAF at Higher Energies, ed. N. 
Isgur and P. Stoler (Newport News, USA, 1994). 

[12] X.-F. Guo, Phys. Rev. D58, 036001 (1998). 

[13] D. Naples, et al, Phys. Rev. Lett. 72, 2341 (1994). 

[14] M. D. Corcoran et al, Phys. Lett. B259 (1991), 209; T. Fields, Nucl. Phys. A544, 565 
(1992). 

20 



[15] G. Curci, W. Furmanski, and R. Petronzio, Nucl. Phys. B175, 27 (1980). 

[16] R.K. Ellis, W. Furmanski, and R. Petronzio, Nucl. Phys. B207, 1 (1982); B212, 29 
(1983). 



[17] J.-W. Qiu, Phys. Rev. D42, 30 (1990); and X.-F. Guo and J.-W. Qiu, [hep-ph/ 98 10515 



[18] EMC Collaboration, J.J. Aubert, et ai, Phys. Lett. B123, 275 (1983). 

[19] A. Bodek, et al, Phys. Rev. Lett. 50, 1431 (1983); 51, 534 (1983). 

[20] NMC Collaboration, P. Amaudruz, et ai, Nucl. Phys. B441, 3 (1995); NMC Collabo- 
ration, M. Arneodo, et al, Nucl. Phys. B441, 12 (1995); Nucl. Phys. B481, 3 (1996); 
Nucl. Phys. B481, 23 (1996). 

[21] E665 Collaboration, M.R. Adams, et al, Phys. Rev. Lett. 68, 3266 (1992); Z. Phys. 
C67, 403(1995). 

[22] K.J. Eskola, V.J. Kolhinen, and C.A. Salgado Eur. Phys. J. C9, 61 (1999). K.J. Eskola, 
V.J. Kolhinen, and PV. Ruuskanen, Nucl. Phys. B535, 351 (1998). 

[23] For example, see F.E. Close, J.-W. Qiu and D.G. Roberts, Phys. Rev. D40, 2820 (1989). 

[24] P. Chiappetta, M. Greco, J. P. Guillet, S. Rolli, and M. Werlen, Nucl. Phys. B412, 3 
(1994); and references therein. 

[25] R. Baier, J. Engles and B. Petersson, Z. Phys. C2, 265(1979). 

[26] J.J. Aubert, et ai, Phys. Lett. B160, 417(1985); M. Arneodo, et ai. Nucl. Phys. B321, 
541(1989). 

[27] H. L. Lai, et ai, Phys. Rev. D51, 4763 (1995); and references therein. 
[28] P. Bordalo, et al, Phys. Lett. B193, 373(1987). 
[29] D. M. Aide, et al, Phys. Rev. Lett. 66, 2285(1991) 



21 



[30] R.J. Fries, B. Miiller, A. Schafer, and E. Stein, |hep-ph/9907567 . 



22 



FIGURES 

FIG. 1. The ol{A,xb) defined in Eq. (g) as a function of xb at (a) Q 2 = 4 GeV 2 and (b) 
Q 2 = 9 GeV 2 . 



FIG. 2. (a) Diagram representing L^ u ; (b) Diagram representing W^ v . 

FIG. 3. Lowest order double scattering contribution to the broadening of the parton c: (a) 
symmetric diagram; (b) and (c): interference diagrams. 

FIG. 4. Lowest order double scattering diagrams which do not contribute to the broadening of 
the parton c. 

FIG. 5. D + (2,z m i n ), -D~(2, z m i n ), and D°(2, z m \ n ) as functions of z m i n . 

FIG. 6. Transverse momentum broadening of pions, A(£j,) 1/ / 3 /A 1 / 3 , as a function of z m i n at 
Q 2 = 4 GeV 2 and (a) x B = 0.2, and (b) x B = 0.3. 

FIG. 7. Transverse momentum broadening of pions, A(£ 2 o) 1 / 3 /^4 1 / 3 , as a function of z m ; n at 
Q 2 = 9 GeV 2 and (a) x B = 0.2, and (b) x B = 0.3. 

FIG. 8. Transverse momentum broadening of pions, 

A(Z 2 )^ /3 /yl 1 /3 ; at q2 = 4 GeV 2 and 

xb = 0.2 with different ^.^^(z). The solid lines are obtained by using the fragmentation functions 
of Ref . |25| , and the dashed lines are obtained by using the fragmentation functions of Ref . 128] . 
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